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Abstract 

We show that two-loop corrections to the finite temperature effective potential 
in the MSSM can have a dramatic effect on the strength of the electroweak phase 
transition, making it more strongly first order. The change in the order parameter 
v/Tc can be as large as 75% of the one-loop daisy improved result. This effect can be 
decisive to widen the region in parameter space where erasure of the created baryons 
by sphaleron processes after the transition is suppressed and hence, where electroweak 
baryogenesis might be successful. We find an allowed region with tan/3 < 4.5 and 
a Higgs boson with standard couplings and mass below 80 GeV within the reach of 
LEP II. 
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1 Introduction 



In recent years a considerable amount of work has been devoted to the study of the elec- 
troweak phase transition in the early Universe. A precise knowledge of it is crucial to 
determine if the very appealing idea of electroweak baryogenesis can be realized successfully 
at all (see for review and references). It is by now clear that the minimal Standard Model 
(SM) has serious problems to generate a sufficient amount of baryon asymmetry at the elec- 
troweak phase transition. The needed CP violation is far too small and in addition the phase 
transition is at most weakly first order for realistic values of the Higgs boson mass. That is, 
the ratio f /T, of the (temperature dependent) vacuum expectation value (vev) for the scalar 
field and the temperature, is significantly smaller than 1 at the moment of the transition. 
Then, erasure of the created B + L asymmetry by unsuppressed sphaleron transitions can 
not be prevented. 

Nevertheless, electroweak baryogenesis may still be successful in alternative models. A 
particularly well motivated extension of the SM is the Minimal Supersymmetric Standard 
Model (MSSM) [Q]. It admits additional sources of CP violation and the presence of 
many extra particles in its spectrum can infiuence significantly the details of the electroweak 
phase transition. A careful study of that transition in this model (see ^ for a study of the 
temperature dependent one-loop daisy improved potential and for previous studies) has 
shown that (if all experimental constraints are appropriately taken into account) a strong 
enough first order phase transition {v/T > 1) can only take place in a very small region of 
parameter space. That region corresponds to large pseudoscalar mass, rriA ^ rriz, small 
tan/3 (~ 2 — 3), third generation squarks as light as otherwise allowed and a SM-like Higgs 
boson with mass just above the present experimental limit {rrih <; 64GeV). 

The purpose of this letter is to go beyond the one-loop daisy improved approximation for 
the finite temperature effective potential in the MSSM in order to see if the small window 
left for electroweak baryogenesis can be significantly widened. Similar studies, where the 
two-loop resummed potential is calculated, have been performed for the SM case and 
sizeable corrections to different phase transition parameters were found. This in turn raised 
some doubts about the validity of perturbation theory and several alternative approaches 
were tried p|. Focusing on 7^ ~ v{Tc)/Tc (T^ denotes the critical temperature, defined here 
by the coexistence of two degenerate minima in the potential), the quantity of interest for us, 
a careful comparison of two-loop resummed perturbative results and lattice simulations (see 
e.g. [0) shows excellent agreement for moderate Higgs masses, indicating that perturbation 
theory is in better shape than previously thought. The inclusion of two- loop corrections 
tends to increase the value of v(Tc)/Tc thus improving the prospects for preserving the 
baryon asymmetry against sphaleron erasure (although the improvement is not sufficient in 
the SM). 

What can then be expected of a two-loop improvement for the MSSM? If the two-loop 
corrections turn out to be small, perturbation theory is trustworthy but the window for 
baryogenesis will remain marginal. If on the other hand very large corrections are found, 
perturbation theory will be suspect, and other techniques should be used. We hope to 
convince the reader that this is a naive expectation. We will show that large corrections 



1 



indeed appear, making the phase transition much stronger, and nevertheless perturbation 
theory could be still considered reliable and under control. 

Before embarking on all the details, let us consider briefly what are the relevant expansion 
parameters in the problem and roughly estimate the size of the corrections expected. In the 
Standard Model the order of the phase transition is determined by transverse gauge bosons 
through its contribution to the cubic term in the effective potential (see [|l| for details). The 
expansion parameter (in the minimum) can be identified as A/ (7^, where A is the quartic 
Higgs coupling and g the SU(2) gauge coupling. The requirement X/g^ < 1 then puts an 
upper limit to the Higgs mass region where perturbation theory can be trusted: roughly 
ml < m^. 

In the MSSM there are two main potential sources of differences with respect to the 
Standard Model: i) the presence of two Higgs doublets makes the effective potential a 
function of two background fieldsQ and ii) contributions of supersymmetric particles to the 
effective potential can have an important effect on the transition. Concerning i), the two 
Higgs doublet potential is in fact tightly constrained by supersymmetry and it turns out 
that the best option for a strong transition corresponds to having only one light Higgs 
doublet (the case with a large pseudoscalar mass). Regarding the new contributions from 
supersymmetric particles it is easily realized that the key role is played by stops. If they 
are relatively light (of course, if all supersymmetric particles are heavy the Standard Model 
result for the corresponding Higgs mass is recovered) they drive the transition: being bosons, 
with a large number of degrees of freedom and sizeable coupling {ht) to the Higgs field, they 
dominate the cubic term in the potential even if their thermal masses are of order ggT. 

In this situation the expansion parameter is \/h1 |Tl|], where ht is the top Yukawa cou- 
pling. Therefore the scale to determine whether the Higgs mass is too heavy to trust pertur- 
bation theory is set now by the top mass Mt while in the Standard Model it was set by My/. 
This is just a naive estimate because the effect of the gsT screening can not be neglected 
and the coefficient in front of X/h^ that would give the proper expansion factor is unknown, 
but it suggests that perturbation theory can now be better behaved than it was in the SM 
for the same Higgs mass. 

However, as we shall see, the situation in the MSSM is radically different from the SM 
case. It is not just a question of whether different effects will make the numerical coefficient 
in front of X/h"^ large or small. The fact that squarks are colored scalars changes the picture 
completely: at two loops there are very important QCD corrections that affect v{Tc)/Tc 
considerably and make the strong coupling gs to appear directly in the game. 

In the next section we review the one-loop daisy improved approximation for the MSSM 
potential This will fix the notation and the interesting parameter range setting the 

stage for inclusion of two-loop dominant corrections, which will be the subject of section 3. 
Results and conclusion will be presented in sections 4 and 5 respectively. Some technical 
details are relegated to the appendices. 

^ Even barring the possibility that other MSSM scalars take non zero vevs at high temperatures one 
should allow for the possibility of a non-zero relative phase between the two Higgs vevs. This case was 
studied in refs. |p^ . The region of parameter space where it can be realized is disconnected with the one 
which will interest us here. 
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2 One-loop resummed potential 



We assume from the beginning that the pseudoscalar mass tua is much larger than Mz- In 
such case only one (combination) of the two Higgs doublets present in the MSSM is light. 
Then, the (light) physical Higgs boson has couplings to vector bosons and fermions of SM 
strength, and in first approximation the experimental limit rrih > 64 GeV applies. This 
simplifies the study of the potential which is then a function of just one background field 
as in the SM. From supersymmetric particles only squarks of the third generation will be 
considered to be at the electroweak scale. This corresponds to the best situation regarding 
the strength of the electroweak phase transition (see ^ for more details). Our model 
reduces then to the SM plus third generation squarks and a quartic Higgs coupling fixed by 
supersymmetry as shown below. 

Working in 't Hooft-Landau gauge and in MS-scheme0, we can write the finite tempera- 
ture effective potential in the form 



where 
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1) The first contribution, eq. (H), is just the tree level part. Note that the quartic Higgs 
coupling is fixed to be 

X = l{g' + g") cos' 2p. (6) 

Here tan f3 is the ratio of the zero temperature vacuum expectation values of the two super- 
symmetric Higgses, {H2) / {Hi) = V2/V1, where H2 is the doublet responsible for the masses 
of up-type quarks and f ^ + f | = (246 GeV)'^. 

2) The second term, eq. (|^), is the one- loop T = contribution. Q is the renormalization 
scale, where we choose for definiteness Q' = m|. The rij's are the number of degrees of 
freedom of the different species of particles, taken negative for fermions: 



nt 



-12 , rii^ =ni^=Q, nw = Q , nz 



rih 



(7) 



and Ci are constants equal to 5/6 for vector bosons, and 3/2 for scalars and fermions. The 
dominant piece of Vi((y9,0) comes from top and stop contributions. 

^In refs. 0, |^ DR-scheme was used instead. We have checked that the numerical difference is irrelevant. 



3 



Next, mi{ip) is the field-dependent mass of the i*^ particle. The field-dependent top and 
bottom masses are given by 

m?(y,) = 1/^2^,2 sin2^ ^ ml{ip) = hilip^cos^p. (8) 

This implies the relations ht^sM = ht sin P and hb^sM = hb cos p. 
The entries of the field-dependent stop mass matrix are 

^IM = ^i + ^tiv) + Dl{v), (9) 
mli^) = ml + m',{^) + Dl{^), (10) 

while the off-diagonal mixing between left and right-handed stops is given by 

m?^^(y?) = -^{AtsmP + fiCOsP)ip = -^Xtif. (11) 

Here rriQ, irtu and At are soft supersymmetry-breaking mass parameters^, is a superpo- 
tential Higgs mass term, and 



1/1 2 

4 V2 
1 /2 



^i(^) = l[l^-^^^^^0w]{9' + 9'^)^'cos2p, (12) 



Dli^) = -(-sin2^^)(/ + ^'Vcos2/3, (13) 

are the D-term contributions to the stop mass matrix. The field-dependent stop masses, 
, are the eigenvalues of the matrix described. 

11,2' 

For sbottoms one has 

^Uv) = ml + m'M + Dli^), (14) 

= ^l + ml{^) + Dl^{^), (15) 

while the off-diagonal left-right mixing is now 



"^Li,'^'^) = -^{^b cos (3 + fi sin (3) if = -^Xbip. (16) 
Here m£),Ab are soft supersymmetry-breaking masses, and 



Dliip) = -i^-- + -sinHw)i9' + 9'Vcos2p, (17) 

Dli^) = l(-lsin2^H^)(/ + ^'Vcos2/3. (18) 
The field-dependent masses for the gauge bosons are given by 

mli^) = \9V, rnli^) = ^{g^ + g'^)^\ (19) 



^For constraints on these parameters we refer the reader to the discussion in 
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and for Higgs and Goldstone scalars by 

ml = 3A(y9^ — /i^ , m^o = m^± = Ay^^ — jj,"^. (20) 

3) The third term in the potential, eq. (^), is the one- loop contribution due to temperature 
effects. Here Jj = J+(J_) if the i^^ particle is a boson (fermion), is the free energy of an 
ideal gas of particles of mass mi{{p), 

J±{y^) = J±{m'/T^) = J^dxx^ log (^1 T e~v/^) . (21) 

4) The last term, eq. (^, is a correction coming from the resummation of the leading 
infrared-dominated higher-loop contributions, associated with the so-called daisy diagrams. 
The sum runs over bosons only. The effect of this contribution is to replace the pure cubic 
term m^{ip) that arises from (^) by (</?). The masses fn^{(f, T) are obtained from the mf{(f) 
by addingfjthe leading T-dependent self-energy contributions, which are proportional to T^. 
We recall that, in the gauge boson sector, only the longitudinal components {Wl, ZL,'yL) 
receive such contributions^ Also note 

= 2 , nz^=n^^ = l. (22) 

The (leading order) thermal polarizations for different particles can be found in ref. 0] and 
are reproduced in appendix A together with some extra polarizations needed for two-loop 
corrections to the potential. We assume here that only SM particles and third genera- 
tion squarks contribute to these thermal masses, the rest being heavy and decoupled. As 
explained in this is the best case for a strong first order transition and in practice corre- 
sponds to the case of heavy gluinos, so that the results found can be interpreted as the best 
guess for a realistic supersymmetric spectrum. 

The contribution of stops to AVdaisy is responsible for the enhancement of7l = v(Tc)/Tc 
with respect to the SM result. The final region in parameter space where 7^ > 1 is determined 
by the interplay between different opposite effects: soft and thermal screening masses tend 
to decrease TZ by screening the pure cubic behaviour of m^{(p) while a large Yukawa coupling 
(thus a large top mass) tend to increase it. For a given value of ht (and tan/3) there is 
a minimum allowed value of mq [the constraint of Ap on the mass splitting of the {t,b)L 
doublet gets weaker for larger mg] and this will correspond to the best case. To determine 
which of the competing effects is stronger one has to study numerically the potential as was 
done in 0, The main conclusion of those papers is that in some regions of parameter 
space the cubic term from stops is effective and actually dominates the electroweak phase 
transition. 

Regions of 7^ > 1 can be found for 135 GeV < < 170 GeV. The best case always 
corresponds to negligible mixing in the stop mass matrix (^f^^ ~ 0), low tan/5 (~ 2 — 3), 

^In general, will be the eigenvalues of the corresponding mass matrix after addition of thermal 
contributions. 

^Longitudinal photons do contribute to the effective potential due to the effect described in the previous 
footnote. 
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a Higgs mass close to the experimental limit {rrih < 70 GeV) and squarks of the third 
generation as light as allowed by the Ap constraint. 

To fix ideas, for rrit = 150 GeV (which corresponds to a physical pole mass Mt = 
156 GeV), rriQ = 70 GeV [this gives Ap{t,b) + Ap{iL,bL) < 0.01], mu = m^-^^ = and 
tan/5 = 2.5, v(Tc)/Tc = 1.02 and so this point in parameter space corresponds roughly to 
the border of the 71 > 1 area. The Higgs mass is ~ 65 GeV. 



3 Two-loop resummed potential 

The complete two-loop resummed potential for the contributions of SM particles can be 
found in refs. 0]. Regarding its effect on the quantity 71 = v{Tc)/Tc, the most important 
corrections are those coming from bosonic setting sun diagrams, which cause logarithms of 
masses to appear. The importance of these contributions was first pointed out by Bagnasco 
and Dine in 0. In particular, the most important effect is due to logarithms of transverse 
vector boson masses, unscreened at leading order. It is enough for our purpose to keep only 
logarithmic SM terms. We will also set g' = for all two-loop corrections. This will simplify 
the analytic results being numerically an excellent approximation. In addition we will use 
throughout a high temperature expansion keeping contributions to the effective potential 
up to order g'^ (where for power counting g represents gauge, Yukawa couphngs or \/X). In 
our formulas, M will represent the (field-dependent) weak vector boson mass and Mi the 
longitudinal, Debye-screened mass (see appendix A): 

M2 = i^V, Ml = M^ + ^g^T^ (23) 



Then, the dominant two-loop resummed SM piece of the potential can be written as 
2 '2/3, Ml 51, M\ 3,^. ^^^^^^ M + 2Ml 



1/(2) _ 9 rj.2 



- log -i^ - - log — + - (M^ - AMI) log 



4°T 8"T/2' 3T 



(24) 



Some comments are in order. Extra logarithmic contributions involving scalar (Higgs) masses 
have been systematically neglected (although these were also kept for the numerical analysis). 
That is a good approximation because near the critical temperature their effective mass is 
small. We have kept logarithms of longitudinal vector boson masses [note that only the 
log(M/T) was kept in Bagnasco and Dine paper]. They are not important numerically 
because the screening of gauge bosons is very effective here but we keep them for consistency: 
similar logarithms will appear for the supersymmetric contribution. Terms like the last one 
in ( p4|) have a linear dependence on <f for small field. As is well known, this linear dependence 
cancels when all resummed two-loop contributions are kept. We shall then add the terms 
needed to insure this cancellation here. These are [0] 

^^sl = ^,MMlT\ (25) 
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There are also non logarithmic contributions involving Qs and ht that can be sizeable in 
principle. They are 0: 



2 r 



5'K 



Here, the constant cb is given by 



lQgi[-\og2--~CB 



^KsM [^^og2 - Cb 



Cb = log(47r) - 



(26) 



(27) 



where je is Euler's constant. This correction affects mainly the numerical value of the 
critical temperature while its effect on TZ is very small. 
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Figure 1: Two- loop graphs involving third generation squarks. 



The extra supersymmetric diagrams^] we have to consider are depicted in figure 1. Note 
that, with our assumptions on the supersymmetric spectrum, R-parity conservation implies 
that squarks appear always in closed loops. Furthermore, there are no fermionic diagrams 
because all R = — 1 fermions are assumed heavy (two- loop diagrams with e.g. a top quark 
and a stop involve also higgsinos or gluinos to close the fermion loop). We label the diagrams 
in a similar way as in refs. 0. By qi we represent a squark of a given flavour q and chirality 
i = L,R. Then, q'j would stand for a squark of different flavour and different chirality. 

The dominant logarithmic contribution (plus linear terms) can then be written, in the 
high temperature expansion (see appendix B for the integral expression), as the sum of the 
following pieces {N^ = 3 is the number of colours): 
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^For MSSM Feynman rules see e.g. |lj 
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In this last contribution we are neglecting stop mixing {fni^^ ~ 0). As was shown in ^ 
for the one-loop potential, this is the best case to maximize TZ. The two-loop resummed 
potential for non-zero stop mixing can be found in appendix C. 

The {z') diagrams do not contribute logarithmic terms but they are needed to cancel the 
linear dependence present in the previous pieces: 

V 

V{z') 



3271' 



;N,T'M{mi^+mr). 



(31) 



From the rest of non logarithmic terms we will keep those involving g"^ and h'f only: 



6V( 
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(A^c + l)(m^-^ + mi^) + h^mi^{mi^ + mi J 



(32) 



where C2 — 3.3025. Here IlgL is the Debye mass for longitudinal gluons as given in ap- 
pendix A. We have neglected throughout the bottom Yukawa coupling, which is small for 
small tan/5. Then terms dependent on m^^^ are not shown. 
The complete terms Sc^V add up to give 



cb 
167r2 



+ m^nf + 3m2 



(33) 



As expected, this term plus a similar term from Standard Model polarizations, combines 
with the one-loop unresummed scalar contribution 

Cb 



327r' 



(34) 
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to give the same result but with rrii —>■ mi (a welcome check of the calculation). In (|33| ) 
]jSUSY jg ^YiQ full thermal mass for squarks while for the Higgs modes it includes only the 
contribution of supersymmetric particles (see appendix A): 

^susY ^ ^susY ^ lh1sm^(3T\ (35) 
2 




^ (GeV) 



Figure 2: Scaled potential at the critical temperature for different approximations: one-loop resummed 
potential (thin solid line) ; potential with two- loop resummed contributions from SM particles only (dashed) ; 
with two-loop resummed supersymmetric non-QCD contributions added (dotted); with all dominant two-loop 
resummed contributions (thick solid). [Mt — 156 GeV , mq = 70 GeV, mjj — rai^^ = and tan/3 = 2.5]. 



4 Results 

In figure 2, different approximations to the scaled effective potential \y{ip)/T^] are plotted 
at the corresponding critical temperature. Parameters are fixed as in the example at the end 
of section 2: top pole mass Mt = 156 GeV, mq = 70 GeV, mu = mi^^ = and tan P = 2.5. 
The thin solid line corresponds to the one- loop resummed potential for which TZ = 1.02. 
The effect of including two-loop corrections from Standard Model particles is shown by the 
dashed line: the transition becomes stronger and TZ increases up to 1.20. This effect can be 
traced back to the presence of corrections of the form —ip'^logip coming from transverse 
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gauge boson modes [see eq. p^)]. A contribution of this form has a direct effect shifting the 
value of the vev at the critical temperature. The minus sign in front implies that the shift 
is towards larger values. To examine two-loop corrections from supersymmetric particles it 
is convenient to split them in two categories: QCD contributions and the rest. Looking first 
at non QCD terms, the dominant effect [see eqs. (^8|-^)] is of the form +93^ log ^ (where 
Tp represents some screened masses). The net positive sign indicates that these corrections 
contribute to make the phase transition weaker. This effect can be seen in figure 2 where the 
dotted line includes non QCD squark contributions. The corresponding TZ falls now down to 
0.97, slightly below the one-loop resummed result. Then, disregarding QCD corrections, the 
convergence of the perturbative series for TZ seems to be much better in the supersymmetric 
case, when the phase transition is driven by squarks, than in the Standard Model (for the 
same Higgs mass). 

The inclusion of QCD corrections changes completely this picture: a look at (pHf ) reveals 
that these corrections are of the form —ip^ ^ogTp and sizeable, implying that TZ will be shifted 
to larger values. The increase in TZ can be as dramatic as shown in figure 2. The thick solid 
line gives the scaled potential with all dominant corrections included. The phase transition 
in this approximation is much stronger and TZ = 1.75 (a correction of 75% over the one- loop 
result!). 

This huge effect makes perturbative validity dubious. Nevertheless, one should keep in 
mind that the effect is very large because QCD corrections appear only at two loops (the 
same happens in the Standard Model but there, the effect is associated with a fermionic loop 
and so the final effect on TZ is much less important) and not at one-loop[|. Then, a large 
ratio of the two-loop to the one-loop result does not necessarily imply that perturbation 
theory is not applicable. To settle definitively this issue one should estimate the size of 
three-loop corrections. However, that is beyond the aim of this paper. Note also that lattice 
simulations would face here the problem that squarks carry color so that SU{'i)c can no 
longer be ignored as in the SM once fermions were integrated out. Lacking an estimate 
of higher order corrections, certainly this two-loop result indicates that the transition is 
probably much stronger than what the one-loop resummed approximation indicated. 

In the previous numerical example we have chosen the parameters to maximize the effect 
on TZ. We explore now how TZ changes when we move in parameter space. Figure 3 is a 
contour plot of TZ (solid lines) in the plane (m(/,tan/5) for rrit = 150 GeV {Mt = 156 GeV) 
and rriQ = 70 GeV. The dashed lines give the Higgs boson mas^ rrih in GeV. Thick lines 
are used to single out the lines TZ = 1 and mh = 64 GeV. The region where electroweak 
baryogenesis could in principle be viable and the LEP I bound on rrih evaded is now sizeable. 
Values of tan (3 larger than in the one-loop case are allowed and Higgs masses up to ~ 80 GeV 
can be reached. The decrease of TZ with increasing mu is just the effect of a growing screening 
in m^^ . For large values of mu the high temperature expansion breaks down and one should 

^This situation is by no means unique and appears in many other cases. Similar examples can be found 
already in the MSSM Higgs sector, e.g. when considering radiative corrections to the mass of the lightest 
scalar Q or QCD radiative corrections to hadronic Higgs decays ]l5{ . 

^To compute the Higgs mass accurately is important to include loop corrections which can be sizeable 
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evaluate two-loop contributions numerically. We cut the x-axis for mu — 140 GeV where 
the high T expansion is still safe. 

Finally, figure 4 shows contour lines of constant TZ and nih in the plane {Mt, tan f3). 
Thick solid lines for TZ = 1 and rrih = 64 GeV delimit a " baryogenesis hole" where the 
electroweak phase transition is strong and the lightest Higgs mass is above the LEP I limit. 
The region corresponds to 125 GeV < Mt < 185 GeV and tan/3 < 4.5. In this plot we have 
set mu — and uiq has been fixed at its lowest value such that the constraint on Ap is 
evaded (fixing for simplicity tan/9 = 2 to compute this minimum of mg). The range of mq 
then varies from zero for Mt = 120 GeV to order 1 TeV for Mt = 180 GeV. The dependence 
of TZ contours on M^ is then easy to understand: for low values of M^ the efi'ect of squarks 
is smaller because the Yukawa coupling is smaller. In this region, increasing Mt one gets 
larger values of TZ. For larger Mt, two effects collaborate to change this tendency. Large 
Mt requires large uiq to avoid problems with Ap so that left-handed squarks are required 
to be heavy. In turn, heavy squarks give a sizeable radiative correction to the Higgs mass 
making it heavier and correspondingly decreasing the strength of the transition. This is 
purely a zero temperature effect. For the finite temperature corrections, heavy left-handed 
squarks are in fact decoupled and do not contribute: TZ in this region has been calculated 
with contributions of right handed squarks only. Therefore, the effect is reduced and TZ 
drops. The net effect is that one ends up with a closed region in the {Mt, tan /9) plane where 
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electroweak baryogenesis may be successful. It is interesting to note that this region requires 
large Mf in a range compatible with the CDF-DO range. Furthermore nih is bounded to be 
below ~ 80 GeV so that LEP II will scan this "baryogenesis hole" in the near future testing 
then one of the crucial ingredients for the viability of electroweak baryogenesis in the MSSM. 



5 Conclusion 

In the MSSM light stops can dominate completely the electroweak phase transition in the 
early Universe. We have shown that the strength of this transition is dramatically affected 
by QCD corrections that first appear at two-loop order. They make the transition strongly 
first order in some non marginal region of parameter space compatible with experimental 
constraints. This is in sharp contrast with the situation in the Standard Model where 
transverse gauge bosons drive the transition which is weak for realistic Higgs masses. 

The allowed region found corresponds to tan /3 small within a reasonable range (tan /? < 
4.5), squarks as light as allowed by electroweak precision constraints and only one light Higgs 
boson, with Standard couplings and mass as large as 80 GeV. Such Higgs scalar is within 
the reach of LEP II so that these region will be explored in the forthcoming years. 

The sizeable increase in the strength of the electroweak phase transition can be decisive 
to open the window for electroweak baryogenesis in the MSSM and should encourage further 
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studies of many different aspects of tlie problem |]T^ still not well understood and necessary 
towards a complete mechanism of electroweak (supersymmetric) baryogenesis. 



Note Added 

In a very recent paper the electroweak phase transition in the MSSM is studied in 
the previously unexplored case mfj < which is shown to be phenomenologically allowed. In 
such region of parameter space the screening of m^^ is then reduced and TZ can be larger than 
1 even at one-loop. Then the allowed areas found in the present paper can be continued 
into the < region improving further the prospects for supersymmetric electroweak 
baryogenesis. 
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Appendix A 

We reproduce here the leading parts of the T-dependent self energies for the different species 
of bosonic particles. Most of these results were presented in ref . . Some errata are corrected 
now (see Ref.[|l^] for details). It is assumed that niA is large so that only one light Higgs 
doublet is present at the electroweak scale. The second doublet is decoupled and does not 
contribute to the thermal polarizations written below. As in ref. terms outside curly 
brackets come from Standard Model contributions and third generation squarks. The rest 
of the supersymmetric contributions is inside the brackets and is written for completeness. 

Note that, at leading order, these thermal masses are equal for particles inside the same 
SU{2) multiplet. They are given by 




n =-{g^ + g' 2) cos^ 2I3T^ + -g^T^ + -g' ^T^ + -hi sin^ /3T' 
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1 9^" 27^ 6 * J ' 

ha g^s ^08^ ^ ' 54^ J 

We complete these formulae giving the polarization for longitudinal gluons 

Ii,, = lglT' + {^-^glT']. 
This last result is needed for a two-loop calculation of the effective potential. 



Appendix B 



We write here the resummed supersymmetric two-loop corrections to the finite temperature 
effective potential in the MSSM framework described in the text. In particular we consider 
g' = hf, = so that in the formulas below we make no distinction between mz and m^y which 
are called M. We assume also that left-right mixing in the stop sector is negligible. Our 
notation is that of Arnold and Espinosa in and the labelling of different contributions 
corresponds to our figure 1. As usual, Nc = 3 counts the number of colours. 



^(4,2) 



9' 



T>ssv {mtL 5 ' ^) + '^ssv {m-f,^ , "^^^ , M) + AVssv {mi^ , m^^ , M) 



- ^{N^ - 1) [Vssvimi^,rni^,0) +'Dssvim-,,,rni^,0) 
+ Vssv{miR,ma^ 0) + Vssv{mi^,mi^, 0) 



ht sin 13+ — cos 2/3 H(m/j, mi^^mi^) + — cos 2/3 H(m,j, m^^, m^^' 



+ {h^sm U{mh,mi^,mi^)+ [h^ sin (3 + -g cos2/3j H(m^, m^-^, m^^ 
= --.gliNl - 1) [Vsvim,.^) + Vsv{m-,^,Q) + Vsviniia^ 0) + 0) 



g'N, Vsv {rriir, , M) + Vsv (m^^ , M) 



(4.4) 
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+ i^hl sin^ (3 + i^f^ cos 2/5^ \ps{mi^,mh) + 2Vs{mi^,m^) + Vs{mi^,m^) 
- ^iVc/ cos2/3 Vs{mi^,mh) +2Vs{mi^,rn^) +'Ds{mi^,m^ 



+ -iVc/ii sin /3 p5(mi^, m^) + 2Vs{mi^, my.) + Vs{mi^,m^) 
where 2)5(7711, is simply the resummed version of 



X>s(mi,m2) = /(mi)/(m2). 

In the previous formulae, dimensional regularization (with n — 1 = 3 — 2e) is used to 
evaluate divergent integrals. Poles in 1/e and i^-dependent terms cancel when counterterms 
are included. In addition, counterterms contribute a finite piece to the potential 0. To 
compute the finite piece due to third generation squarks note that the counterterm potential 
can be written now as 

+ i(3 - 2e) {3g\6Zg. + 6Z^.) + g'\5Zg,. + 5Z^.))^ - 6^%,,/;^(5Z,2^^^ + 5Z^.)^ 



+ ^ip\5Zx + 5Z^2), 



(36) 



with 



' ^//3 = -^[l + <^^-log4)]. 



The MS renormalization functions 5Z have no finite part, so the finite (te-independent) 
piece of 5V is given by the (3 — 2e) and JJ-^ terms. The required contributions of third 
generation squarks to the renormalization functions are 



6n„Zn 



18" levr^e' 



;9 



5g',Zl^2 



0. 



(37) 



The third equation follows from the fact that squarks are scalars while the fourth is a conse- 
quence of our assumptions about the supersymmetric spectrum, that has no light R = — 1 
fermions. The finite contribution is then 



6V 



167r2 96 



(3J 
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Appendix C 



To include the effect of the left-right squark mixing appropriate mixing angles, at for stops 
and «(, for sbottoms, should be defined. When temperature corrections are included in 
the mass matrix for the zero Matsubara modes, the corresponding mixing angles will turn 
out to be temperature dependent. This effect has to be taken into account when writing 
the resummed effective potential and zero modes should be treated independently. It is 
then possible to generalize the previous expression for the effective potential writing it in 
terms of unresummed "D's plus extra terms for the resummation of zero modes. A similar 
situation arises with the photon-Z mixing in the Standard Model. In the same way that 
this complication is circumvented in the Standard Model just by setting g' = 0, a. reasonable 
assumption for the squark case would be to approximate the squark thermal masses by its 
universal QCD contribution Ag'^T'^/d and neglect the rest. In this way the thermal correction 
to the squared mass matrices would be diagonal and the mixing angles would be temperature 
independent. In this approximation, the generalization of the previous expression for the 
potential is straightforward. Writing 



ii — cos attL + sin atiji 
t2 — cos attn — sin aft l 



cos abbi + sin q;(,6r 
62 = cos abbfi — sin a^bL 

we obtain (using the abbreviations q = cosa^, C2/3 = cos 2/3, etc) 



(39) 



9 r 

V(a') = [ctl^ssv{'mi^,'mi^,M) + s^Vssv{mi^,mi^, M) + 2cfs^Vssv{mi^,mi^, M) 

V"c [4Vssv(mi^,mi^,M) + stVssv(mi^,mi^, M) + 2clslVssv(mi^,mi^, M) 



2 ' 



+ V 



T^ssv [mi^ , m^^ , 0) + Vgsv {mi^ , , 0) 



'ssv {mi2 > "^£2 > 0) + T^ssv {mi^ , m^^ , 0) 



(p') 



+ 



+ 



9^ , , X 



9 



H(mfe,mt^,m£j 



9 

{hts} + —StC2fs)(p + ht\/2{AtSi3 + i^cp)ctSt 



^ipCbSbC2p + -^{AbCp + IISp)c2b 
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+ 



+ 



^2 

—(pclc2f} - hby/2{AbCf3 + iJ,Sf3)cbSb 

'g2 

-^Vslc2l3 + hbV2{AbC/3 + IJ,Sp)CbSb 



1 2 



Nc 2(9 



ht {Atsp + ficp) B.{m^, mt^,mtj + hbiAbCp + i^sp) B.{m^, mi^,mi^) 



■^{h^sl + ^C2i3)ctCb + htXtStCb - hbXbCtSb 



+ 



+ 



+ 



^{hlsp + yC2/3)CtS6 + htXtStSb + hbXbCtCb 



;^(^t«J + ^C2p)stCb - htXtCtCb - hbXbStSb 



1 2 



-j={hfSf^ + —C2p)stSb - htXtCtSb + hbXbStCb 



where we keep hb only when it appears in combination with Xb- 

y _ 1„2,..2 



'(4,2) 



—ai{N'^ - 1) PsyK^, 0) + 2^5y(m-,^, 0) + PsyK,, 0) + 0) 



+ '^s{mb^,m-J 



Nc 



CbS^Vs{mi^,mi^] 



X g^c23 

- \ ^ [cl[Ds{mh, m^J - Psl^;^, m^J] + sl[Vs{mh, m^J - ^^(m;^, m^J] 



2 



(TVc + 1) Vs {mi^ , m^^ ) + Vg {mi^ , m^^ ) - 2NcC^ Vs{mi^ , mi^ ) 



+ -r^c(2 - Nc) CtCbVs(mt^,mi^) + CtSbVs{mi^,m~^) + StCbVs(mt^,miJ 
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1 



+ ^/iV,(iV, + l) 
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